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Abstract 



Using the analytic, global solution for the rigidly rotating disc of dust 
QQ I as a starting point, an iteration scheme is presented for the calculation 

. of an arbitrary coefficient in the post-Newtonian (PN) approximation of 

this solution. The coefficients were explicitly calculated up to the 12* PN 
, level and are listed in this paper up to the 4"^ PN level. The convergence 

of the series is discussed and the approximation is found to be reliable 
even in highly relativistic cases. Finally, the ergospheres are calculated 
CJ> , at increasing orders of the approximation and for increasingly relativistic 

situations. 

U 
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1 Introduction 

■ The study of a three-dimensional, arbitrarily rotating fluid in general relativity 

5-H ' cannot but rely on the use of numerical computation for its results. An impor- 

^ tant step along the way toward a proper application of numerical methods and 

toward an in-depth understanding of numerical results is the rigorous consider- 
ation of a simpler, but related, problem. Neugebauer and Meinel undertook this 
task in modelling a uniformly rotating disc of dust. The corresponding analytic, 
global solution to Einstein's field equations, which was worked out in the series 
of papers Q, |^ and by utilizing the 'inverse scattering method' known from 
soliton theory, can be expressed using hyperelliptic integrals. 

In this paper, the above mentioned global solution will be used as a starting 
point from which to derive an iteration scheme for the calculation of an arbitrary 
term in the post-Newtonian (PN) expansion of the solution. Such an expansion 
amounts to the analytic analogue of numerical work presented by Bardeen and 
Wagoner Q in a paper that handles the rotating disc in great detail. 

Given that the global solution for the uniformly rotating disc of dust is known, 
the question arises as to why one would consider the problem using the PN 



1 



approximation. The reasons are threefold. First, the complex nature of the 
solution leads to fairly long computing times, particularly if the relativistic 
parameter /i is to be varied. Using the PN approximation, one could speed 
up calculations dramatically. Secondly, one is presented the rare opportunity 
to check the accuracy of numerical work against its analytic analogue, thus 
allowing one better to determine the limitations of these numerical techniques. 
Lastly, since one can obtain an arbitrary term in the PN expansion, one can 
look at the convergence of the series for various applications and, moreover, see 
if the approximation can be used to study relativistic phenomena, where its 
validity is not a priori clear. 



2 The Global Solution 

The form of the metric to be used here is taken directly from Q and reads 

ds^ = e-2C^[e2'=(dp2 + dC^) + p^d^^] - e^^'idt + ad^f (1) 

in Weyl-Lewis-Papapetrou coordinates. [| The metric represents a stationary, 
axially symmetric spacetime and the three metric functions C/, k and a depend 
only on p, C, the relativistic parameter 

002 2 -2Vo 

^^2^ 

and the coordinate radius of the disc po- The other quantities appearing in 
eq. (|^) arc the constant angular velocity and the "surface potential" Vq = 
J7(/3 = 0, C = 0), which is closely related to the redshift at infinity. The relativistic 
parameter /i runs from /i = in the Newtonian limit through to /i = /xq = 
4.62 ... in the extreme Kerr limit. The four-velocity of a particle in the disc has 
only two non-zero components: u* — e~^" and u"^ = ne~^° (both components 
are independent of the radial coordinate p). As mentioned above, the metric 
depends on only two parameters, whence the profile for the surface mass density 
cr{p) cannot be chosen freely, but is instead automatically determined by these 
parameters, i.e. a disc chosen to have the values fJ- — p and po — Po can have 
but one angular momentum J and but one total gravitational mass M and the 
matter will needs be distributed according to a{p) as given by eq. (23) in 
A specific example of a mass density profile will be presented in section [4 .21 . 

The vacuum field equations are equivalent to the complex Ernst equation (see 
[|) 

5i(/)A/ = (V/)2 with / = e2^+*6 (3) 
and b p ~ a,^, 6^^ — a^p. (4) 




^ Units have been chosen in which G and c are equal to one, except in eq. (y), where the 
factor w as included to show explicitly the relation fj, oc for reasons to be discussed in 
section 

^One could just as easily have chosen to consider a disc of mass M' and angular momentum 
J' and would then have been automatically led to the values for fi, po and cr(p) corresponding 
to this situation. 
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The operators A and V in the above equation are three-dimensional. Taking 
into account boundary conditions on the disc, an asymptotically flat solution 
can be found for the Ernst potential /, which suffices to determine all three 
metric functions since a can be found using eq. (Q) and k is related to the 
other two functions via a line integral as well. The potential / is given by the 
expression 

^■'^ Jo 2Xim)Wiim) 



{Xin)-X,)iXin)-X2) 
2X{n)Wi{n) 

H{X - Xi){X - X2) 



dn 



dX (5) 



W 

h+h- h (6) 



with 



W = W^W2, W, = -^{X-C/pof + {plpo)\ (7) 

w, = v/^F^ToTW, i^^=»±^, (8) 



X, ^ -^tzJi „a X, ^ ^-i±t. (9) 

A negative sign appearing before the root of a complex quantity indicates that 
the real part is to be chosen to be negative. What is meant by X{in), X{n), 
Wi (m) and Wi {n) in eq. (0) is that the variable substitution 



icosh(m) _ /— icosh(n) _ 

X^-^^y-1 or X^^^-^-1 (10) 

is to be carried out. The endpoints of integration TOq and can be determined 
from the Jacobi inversion problem 



dm /'"'' dn 



and 



with 



2X{m)Wi{m) Jq 2X{n)Wi{n) 
dm /■"'' dn 



2Wiim) Jo 2Wiin) 



(11) 
(12) 



HdX , r HXdX , , 

and - / (13) 



W J_, W 

and where the path of integration in eq. ( p^ ) , as with the integral in eq. (F 
is along the imaginary axis. The above equations, which can be found in I 
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differ from the ones presented there in three ways: 1) The expression for ln(/) 
was manipulated algebraically, making use of eqs. (|l^), ( p^ ) and (|l^). 2) In 
each of the integrals Ii and I2, the aforementioned variable substitution was 
carried through. 3) The definition for W2 in this paper differs from that of 
by a factor ^. 

In discussing the Newtonian limit of the solution in ||^, Neugebauer and 
Meinel make note of the fact that the solutions can be expanded in a power 
series 

oc 
n=l 

such that the coefficients are elementary functions. The following section presents 
a method for determining these coefhcients. 



3 The Iteration Scheme 

A comparison of eqs. ( ]lO| ) and (H) shows that 

X{m ^0)^X1 and X{n = 0) = X2 (15) 
hold. This suggests the rearrangement of eqs. (0) and (|l^) to form 



- X{m)-X2 , ^ r X{n)-X2 , ^ 
I -, — ^ -, — r dm + / an = v — Xou , (16) 

X(m)~Xi , r" X(n)-Xi , , , 

/ -^dm+ / -i-dn = v - Xiu, (17) 

Jo 2Xim)Wi{m) Jo 2X(n)H^i(n) ' ' ^ ^ 

which serves to "decouple" the original equations. What is meant here by decou- 
pling, is that an iteration scheme for the determination of rUa (or equivalently 
rib) as a power series in fj, does not require the simultaneous consideration of 
two equations in its ultimate step. Eq. ( p6| ) yields up the ultimate term for to^ 
and eq. ( p^ for rib. 

Expanding the integrands in eqs. (|l6|) and ( pT[ ) about the points m — and 
n = respectively, one is left with trivial integrals 



00 00 



"^ai I m^' dm + ^ Ci / n^* = v ~ X2U (18) 

'^d, / m^'dm + ^h, I n^'dn = v - Xiu. (19) 

Making use of the fact that nia and rib are both of the order 0{i.i), 
derive the following iteration formulae: 



1 

= — 

ao 


{Vifi^ + ■ ■ 


•^This follows from eqs. 


Q) and (0 



along with the fact that the coefficients ai, hi, Ci and 
di are all of the order 0{y^^, that u and v are of the order 0{fj,'^) and that Xi and X2 are 
of the order 0(1/^). 
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-ai — 02 — - — 



,2fc-l 



2fc- 1 



-Ci- 



Cfc-l 



2A:- 1 



(20) 



and 



1 



2k\ 



-Oi — 02 — = — 



,2fe-l 



TO 



— • • • — Ofc_i— — 



2fc- 1 



dk-1 



,2fc-l 



2fc- 1 



where TOq^ and are defined by 



ma^ma,+0{p'^+') and = n6, + 



(21) 



(22) 



and Uj as well as Vj are defined below in eq. (psj). 

It proves useful to augment the collection of integrals u and v by defining a 
third integral 

^^'^^ (23) 



Upon introducing the oblate spheroidal coordinates ^ and rj 

P = Pov/(l + e')(l-'7'), C = < e < 00, -1 < 77 < 1 (24) 

one can come up with very simple, closed-form expressions for an arbitrary term 
in the series expansions of u, v and w about the point fi = 



u = ^ Ujfj,^-' , V = ^ Vj^^^ and w — ^ wjfi 



2j 



These expressions are given by 



arccot(^) 



(25) 



(26) 



arccot(^) 





arccot(^) 



tan' {g)P^^-^ dg 



a J -I 



(27) 



(28) 
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with (3 

7 

and a 4 



1 - (1 + e){l - V') sin'ig) ~ tan^g), 



^(2j + !)!!■ 



(29) 
(30) 

(31) 



An expansion of the integrals Ii and I2 yields an expression in terms of ma^. 
and nhj. . ^ These can in turn be expressed as a series in /i by expanding the 
coefficients a^, bi, Ci and di about the point /i = 0. The integral I3, which can 
be written in the form 



(32) 



can easily be converted into a power series in [i by expanding X\ and Xi about 
the point /i = 0. Thus a means of determining the coefficients /„ of eq. ( p^ 
and representing them in terms of Uj, Vj and Wj has been found. 



4 Results 

4.1 The Ernst Potential 

The first eight coefficients in the expansion of the Ernst Potential are 

/i 
h 

h 



and 



1 2 

2^1 -"^1' 



/4 = 

h = 

/6 = 

/7 = 



1 



\/2wiWi H 

\/2 



1 1 3 2 

U2 + 77^1 + -u^ - UxWx - Wi , 

2 

"24" 



- 3wi - 24i;2 + 8^77Ui3 - 24ui2vi 



-12ui wi + 24 w\ 



-2 Ml Vi 

1 



(33) 
(34) 

(35) 
(36) 
(37) 

(38) 



(39) 



*One should be careful to nlinosn the same sign for square roots in Ii and I2 as were chosen 
upon expanding eqs. ( [ll| ) and (|l2[). 
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/s An 



- 24 wi'^ ~3vi 



48" 

+30 Ml vi — 24 ui wi — 4 
48 vi U2 



24 W2 - 8 "1^ + 48 ui^ vx 
\r]^u\ V]^ — 24 U]^ v\ 



(40) 



Because Wj, and are all real quantities, one can see immediately that 
the coefficients /„ are alternately real and imaginary. Since the relativistic 
parameter is proportional to this series clearly exhibits the structure of 

the PN approximation, in which terms occur in pairs, the first of which is 
imaginary and of the order ©(l/c^™^^), m e N, and the second of which is real 
and of the order ©(l/c^™). Note that /i represents the Newtonian limit, i.e. 
the first truly post-Newtonian contribution is given by m = 2. Although the use 
of the expressions Uj, Vj and Wj provides a fairly succinct notation for the /„, 
the coefficients still quickly become unwieldy with increasing n. For example, 
the expression for /24 would fill approximately 30 pages. 

When the full expressions for mi, wi, w\ and W2 are substituted into eqs. (p3|), 
55|), ( |3^ ) and (^, one obtains, using the abbreviation arccot(^) = x, the first 
four coefficients in the expansion of e^^ , e^^ = 1 + /2n-iM", 



/i 

h 



1 r 



127r 



(54e'^X^ + IOStt^^x- 108C^x + 36C^x^ - IOStt^^ - 
54^2 - 36Cx + 6x^ - 55 7rC + 9 7rx)??* + ( - 36^'*x^ 
90 TT ^4 ^ + 72 X + 24 X' + 90 ^ - 72 ^ X - 36 e 
12x' + 42^e-67rx)?7' + 6^"x' + 9^e"x-12C'x 



(41) 



90 7r^2^ 
-24Cx 



6^e'x + 6e' + 12ex + 6x' 



37r^ - 37rx 



(42) 



nO^^ X^ + 9225 7rC'^ ^ - 3465 tt^ ^ - 2430 x^ + 810 x^ 



2 c5 



2430Cx- 4725 7r^ Cx 

3 , 



9225 TT - 15900 vr ^ + 3375 tt 



1807r3 

- 18450 TT X + 1 1025 TT x^ + 3465 tt^ C 

- 1620^3 X^ + 270 X^ 

- 810 ^ + 3570 TT^ - 1575 71^^^^ + 810 x - 270 C X^ + 30 x^ 

- 4875 TT - 2670 vr C X + 135 TT x^ + 693 TT^ ^ - 75 TT^ X + 320 vr) rf 
( - 810 X^ - 11025 vr^^ x^ + 4725 Ti'^fx + 2430^^ x^ 

- 270 X^ + 22050 tt x - 9405 tt x^ " 4725 tt^ _ 2430 x 

- 4725 TT^ ^"^ X - 540 f ^ + 450 x^ - 1 1025 tt + 13560 tt ^ 



- 1215 TT ^2 x^ + 810 - 3150 TT^ 



675 TT^ f X + 270 X - 450 ^ x^ 



90 x"^ 



4155 7rf2 + 690 7r^X + 45 7rx^ 



45 -K^i-^h-K^ X)Tf 
fi ,,3 



+ (270 x^ + 3375 tt(^x^ - 1575 tt^ x - 810 x^ - 450 ^"^ x 



6750 TT^^ X + 1215 TT^ x^ + 1575 + 810^^^ X - 675 tt^ C X 
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- 900 + 90 ^2 + 3375 tt - 1980 tt x - 855 tt _ 270 
- 150 TT^ + 405 n^fx- 450 x - 90 ^ + 90 + 765 tt 
-630 7r^X- 135 7rx^ - 315 tt^ ^ + 45 tt^ x)??^ - 30 ^^x^ 
135 TT ^6 x' + 75 TT^ ^6 X + 90 e x' + 90 ^4 ^3 ^ 270 TT ^5 X + 45 TT 

! c5 _ on t4 ^, _ /( r „2 a4 ,^ _ -1 on t3 ,,2 _ on t2 ^,3 _ 1 or „ t4 



- 75 TT^ - 90 X - 45 TT^^^x- 180 x' - 90 x' - 135 tt 
+ 135 TT ^2 X^ + 30 ^3 + 70 TT^ f - 45 TT^ ^2 X + 90 ^2 X + 90 ^ x^ + 30 X 

- 457r^^- 907r^X- 457rx^ + 157r^$ + 157r^x] (■ 

nd 

/7 = Y^^^[(34020C8^^^ 19089007r^8^3_;^gg245 7r2^8^2_^3g(3g(3^7^3 
+ 45360 £^ - 675675 t:^£^y-- 5726700 tt^^ + 30139: 



,3 

(43) 



- 45360 f x" - 675675 n^fx- 5726700 tt^^ x^ + 3013920 tt x^ 
+ 392490 TT^ X - 326340 tt^ x^ + 204120 x^ - 136080 x^ 
+ 22680 x'^ + 675675 tt^ + 5726700 n^^x- 1261260 tt^ x 

- 7132860 TT^^ x^ + 1544760 TT^^ x^ - 196245 tt^ + 521850 tt^^^ x 

- 136080 X + 136080 x^ - 179550 tt^ x^ - 45360 x^ 
+ 5040 ^2 x"^ + 1036035 tt^ - 1908900 tt - 727650 tt^ x 
+ 5223960 TT^^'x - 2560740 tt^^ x^ + 272160 tt^^ ^3 _^ 34920 

- 195510 TT^ ^4 - 45360 x + 168294 tt^ x + 22680 f x^ 

- 44100 TT^ ^2 _ ^ ^3 _^ ^4 _^ 442355 ,^3 ^3 - 1105020 tt 

- 132300 TT^ X + 1123500 tt x - 237300 tt ^ x^ + 7140 tt x^ 

- 10549 TT^ + 16758 tt^ ^ x - 1365 tt^ x^ + 45477 tt^ ^ - 107520 tt ^ 
+ 35840 TT X - 3675 tt^ x) + ( - 45360 x^ 

- 3013920 TT^® x^ + 326340 tt^ ^2 _|_ i^i^^q ^7 ^3 _^ ^261260 tt^ x 
+ 9041760 TT^^ x^ - 3220560 TT^^ x^ - 652680 tt^ x - 219240 tt^ 

- 272160 x^ + 30240 Sf" x" - 1261260 tt^ + 2134440 tt^ x 



- 9041760 TT X + 7672560 tt x^ - 619920 'k^'^x^ + 326340 tt^ 
+ 181440 X + 397320 tt^ x - 655200 tt^ x^ - 60480 x^ 
+ 13440 ^^x^ +3013920 TT^^ - 17140207r^^^ + 10584007r^ x 

- 5683440 TT X + 954240 n^^x^ + 122640 tt x^ - 178080 tt^ 

- 45360 Sf" + 854952 tt^ x - 227640 tt^ ^^.2 ^ 3Q240 ^2 ^2 

- 13440 ^x^ + 1680 X* - 599172 tt^^^ _^ 1231440 tt^^ - 370160 tt^^^ 
+ 147000 TT^ X - 156240 tt ^ x^ + 1 1760 tt x^ - 249732 tt^ 



i4/uuu7r 5 X ~ J^ooz4U7r^ X i-ii'DU7rx — z^y/ozTr 5 
. 169848 TT^ ^ X - 10500 tt^ x^ - 40908 tt^ ^ + 35840 tt ^ + 35840 tt x 
+ 2100 TT^X- 17920 Tr^)??^ + (22680 ^^x^ + 1544760 tt^^ x^ 

- 179550 TT^ ^8 x^ - 90720 x^ - 30240 x* - 727650 tt^ ^ 

- 4634280 TT^^ x^ + 619920 tt C*^' x""* + 359100 tt^ x + 136080 x^ 
+ 655200 TT^ ^6 -^2 _^ 90720 x^ - 5040 i'^ x^ + 727650 tt^ ^ 
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+ 4634280 TT X - 1058400 tt^ x - 1662360 tt^^x^ - 594720 tt 

- 179550 TT^ ^6 - 1077300 tt^ f"^ x - 90720 x - 90720 x^ 
+ 779940 TT^ ^4 _^ ^QQgQ ^3 ^3 _^ ]^0080 _^ 815850 ^5 

- 1544760 irf- 396900 tt^ x + 1464960 tt x + 884520 tt x^ 

- 196560 TT x^ + 422100 tt^ + 22680 + 30240 x 

- 886620 TT^ X + 146160 tt^ x^ - 5040 x^ - 10080 ^ x^ + 2520 x^ 
+ 189630 TT^ ^3 - 422520 TT ^3 - 25200 tt^ ^ - 289800 tt x 

+ 130200 TT ^ x^ - 2520 tt x^ + 204330 tt^ _ gg^QQ ^2 ^ ^ _ ggg ^2 ^2 

+ 630 TT-"^ e + 630 TT^ x)??'* + ( - 5040 f x^ - 272160 tt^^ x^ 

+ 44100 TT^ x^ + 20160 x^ + 13440 x'' + 132300 tt^ x 

+ 816480 'kCx^ + 122640 tt^^ x^ - 88200 -k^Cx- 227640 tt^ x^ 

- 30240 S,^ x^ - 40320 x^ - 10080 x* - 132300 tt^ 

+ 147000 TT^^^x- 816480 TT^^x- 186480 TT^^ X^ + 196560 tt^^ x^ 
+ 44100 TT^ ^6 + 20160 X + 367080 tt^ x - 146160 tt^ x^ 
+ 40320 X^ + 20160 C^X^ + 272160 TT^^ - 102900 tt'"^ + 5040 tt^"^ x 
+ 25200 TT^ X - 302400 tt x^ - 35280 tt x^ - 139440 tt^ 

- 5040 - 13440 X + 140840 tt^ x - 10080 x^ + 22680 tt^ x^ 
+ 1680 x^ - 2660 TT^ f + 58800 tt + 105840 -k^x- 2520 tt^ x 

+ 352807r^X^ - 11760 ttx^ - 28980 tt^ _ 24350 tt^^ x + 46207r^ X^ 
+ 420 TT-'' C + 420 TT^ x) + 420 X^ + 7140 TT X^ 

- 1365 TT^ f x^ - 1680 C x^ - 1680 x^ - 3675 tt''' x - 21420 tt^^ x^ 

- 11760 TT x^ + 2730 7r2^^x + 2520 C^x^ + 10500 tt^^'^ x^ 

+ 5040 x^ + 2520 x^ + 3675 tt^ + 21420 tt x - 2100 tt^ x 
+ 18900 TT x^ - 2520 tt x^ - 1365 tt^ - 16310 tt^ x - 1680 x 

- 5040 x^ - 630 TT^ x^ - 5040 x^ - 1680 x" + 875 tt^ 

- 7140 7r^5 + 630 TT^ x - 2520 tt x + 16380 tt^^ ^2 ^ nygOTr^^ ^3 

+ 5810 TT^ ^4 + 420 ^4 + 1680 x + 1190 tt^ x - 4620 tt^ x^ 

+ 2520 ^2 x^ + 1680 C X'* + 420 x* - 595 tt^ - 4620 tt - 420 tt'' x 

- 138607r^^x- 13860 TT^x^ -4620 7rx^ + 1155 77^^^ + 2310 TT^^X 

+ 1155 TT^ X^ + 525 7r^^ + 525 77^x1- (44) 

Calculating the appropriate line integrals, one can determine the remaining 
metric functions a and k. Using the expansions 

00 00 
a=^a2„M^2n+i)/2 e2'= = l + ^jr2„-iM", (45) 

n=l n=l 
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one finds the expressions 



02 = 



04 



\/2/90 



06 



1 - 



(15 ex-i5e+isex-rd^+3x)v^ 



-^ex+^e-^ex+c+x , (46) 



l-rf 



(180 X' - 360 X + 252 x' + 180 - 384 ^ + 108 e X 



2 



1-ry^ 



( 1 2960 x^ + 7875 TT - 45045 tt^ x - 38880 x^ 



727r2 

+ 132 ^2 _ 24 ^ ^ + 36 ^2 _ g4^^4 ^ ( _ 72 ^6 ^2 ^ ^^44 ^5 ^ 

+ 216^^X^ + 1357r^^x- 72^^ + 360^^x^ - 1357r^^ + 1627rC^x 
-216^^ - 144^x + 72x^ - 1177r^ + 277rx - 64)77^ 
+36 e - 72 ^5 ^ - 36 x' - 27 TT X + 36 - 36 x' + 27 TT ^ 
-18 7re'x + 36e' + 72^X + 36x' + 97r^ + 97rx-64], 

14407r3 

+ 25920 - 15750 tt x + 15300 tt x^ + 45045 tt^ + 38880 x 

- 97020 TT^ X - 64800 x^ + 17280 x^ + 7875 tt - 25350 tt^^x 
+ 9450 TT x^ - 12960 + 82005 tt^ + 51840 x - 66150 tt^ x 

- 30240 x^ + 4800 x^ + 10050 tt^" - 1 1850 tt x + 1620 tt x^ 

- 12960 e"^ + 42819 tt^ + 1296O ^ _ 14700 tt^ ^ _ 4329 ^ ^2 

+ 480 x^ + 3275 tt - 2250 tt ^ x - 405 tt x^ + 5619 tt^ ^ - 525 tt^ x 
+ 1280 tt) + ( - 12960 x^ - 7425 tt x^ 

+ 51975 TT^ X + 38880 x^ - 8640 ^'^ x^ + 14850 tt x - 8100 tt x^ 

- 51975 TT^ - 38880 x + 94500 TT'^fx + 30240 x^ + 1 1520 x^ 

- 7425 nf + 11250 tt x - 630 tt x^ + 12960 - 77175 tt^ 

- 34560 X + 47250 ^4 X - 15840 X^ + 8640 ^2 X^ - 3150 TT ^4 

- 4530 TT ^3 X + 540 TT x^ + 12960 f - 26145 tt^ + 4320 ^ 

+ 4500 tt'^^^x- 7200 ^ x^ + 1440 x^ + 4335 tt - 930 tt ^ x + 495 tt x^ 

- 225 TT^ ^ - 225 TT^ X - 640 tt) t?* + (4320 x^ 

+ 2025 TT ^® x^ - 14175 tt^ ^® x - 12960 x^ - 2880 x^ - 4050 Tr(Jx 

+ 1980 TT^s x^ + 14175 TT^ + 12960 C*^ x - 18900 tt^ ^« x + 1440 x^ 

- 5760 x^ + 2025 irf- 2610 tt x + 4590 tt C'' x^ + 14175 tt^ 

- 4320 - 3375 tt^ x + 5760 x + 12960 x^ + 2880 
+ 6307r^'' - 990 7rC^x + 5580 7r^2 -^2 _432o^3 _ gQ^2 ^3 

+ 1350 TT^fx- 7200 X - 1440 ^ x^ + 1440 x^ - 3375 tt f 

- 2430 TT ^ X + 945 TT x^ - 720 TT^ ^ - 640 7r)r?^ 

-480 x^ + 405 TTfx^ + 525 tt^ x + 1440 x^ + 960 x^ 

- 810 TT X + 900 TT X^ - 525 tt^ - 1440 x + 300 tt^ x 



(47) 
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- 1440 + 405 TT - 1530 TT ^5 X - 450 TT - 125 TT^ ^5 + 480 

- 225 n'^S.'^x- 1440 x^ - 960 £^ x^ + 630 tt + 90 tt x 

- 540 TT ^2 ^2 _^ 220 + 439 ^3 + 1449 _ 30 tt^ ^2 ^ _^ ^44Q ^ ^2 

+ 480x^ +405 7r^2 _^810 7r^X + 405 7rx^ - 30 7r2 ^ - 30 7r2 x - 640 tt , 
0, 



^ - " (9e'x'-18C'x+10C'x' + 9C2-14^x + x'+4)r?2 



1-r;^ 



1-772 



27r2 

-C^x' + 2e'x-2€^x'-f -2Cx-x'+4 
1 

"18^ 

+ 33^2-6ex + 9x'-16)r/4+ (- ISC'' X' + 36 ^^X + 54^4 x'-lSe^ 
+ 90C'x' - 54^2 - 36CX+ 18x' - 16)??' 

+9C'x' - IS^'x - 9C^x' +9^^ - 9C'x' + 9^' + 18Cx + 9x' - 16 



(45 e X' - 90 C^ X + 63 e X' + 45 - 96 ^ + 27 ^2 x2 



- [1 - r?'] [ (7290 e X^ - 7245 ^ f x' - 29160 x^ + 16200 ^ 

+ 14490 TT^ C^ X - 16020 tt^ ^'^ x^ + 43740 x^ - 55080 x^ 

+ 10620 ^'^ x** - 7245 ir"^ + 27210 7r2 ^ - 29160 x - 10350 tt^ C* x^ 

+ 68040 x^ - 28440 x^ + 1800 ^2 / + 7290 C'' - 1 1 190 7r2 

+ 14934 TT^fx- 35640 x - 1620 tt^ f x^ + 26460 ^2 ;^2 _ 2520 ^ x^ 

+ 90 x'' + 6480 ^2 - 5389 n'^ + 1734 tt^ ^ x - 10080 ^ X - 45 tt^ x^ 

+ 720 x^ - 256 TT^ + 1440) r;^ + ( - 8910 f + 8775 tt^ x^ 

+ 35640 C^ x^ - 21240 f x* - 17550 tt^ x + 22500 tt^ x^ 

- 53460 x^ + 64440 x^ - 16020 x"^ + 8775 n"^ - 29070 7r2 ^ 
+ 35640 C^ X + 17730 x^ - 65880 x^ + 33480 f x^ - 3960 x" 

- 8910 C* + 6570 7r2 - 13890 tt^^\ + 23400 x + 4500 7r2 ^2 ^2 

- 17460 ^2 x2 + 4680 ^ X^ - 270 x^ - 720 ^2 + 495 - 930 7r2 ^ x 

- 1440 C X + 495 7r2 x^ - 720 x^ - 896 tt^ + 1440) r?^ 

+ (1710 f x^ - 1575 TT^ C^ x^ - 6840 (J x^ + 5400 x* + 3150 tt^ x 
+ 10260 x^ - 5580 tt^ x^ - 9720 x^ + 5940 x* - 1575 tt^ 
+ 3390 7r2 ^5 X - 6840 i^x- 3240 S,'^ x^ - 2970 tt^ x^ - 4680 x^ 
+ 2520 ^2 / + 1710 ^4 _^ 2190 -k^ + 14040 x + 1410 7r2 ^ 
+ 1620 7r2 ^2 ^2 _ ggQQ ^2 ^2 _ -^gQQ ^ ^3 ^ 270 x^ - 6480 ^ ~ 855 
+ 10080 e X - 270 7r2 C X + 585 7r2 x^ - 720 x^ - 1440 - 896 n^) if 
-90 x" + 45 7r2 ^8 x^ + 360 x^ - 360 x^ - 90 tt^ C x 

- 540 S,^ x^ + 540 TT^ x^ + 360 x^ - 540 x* + 45 tt^ - 570 tt^ x 
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+ 360 X + 1080 f " - 90 TT^ ^2 _ ggQ ^3 ^3 _ ggg ^2 ^4 _ 9Q ^4 

+ 30 TT^ ^4 - 1800 X - 150 TT^ ^3 X - 180 tt^ ^2 ^ qqq ^2 ^2 „ ggg ^ ^ 

- 90 x'' + 720 $2 + 405 tt^ ^2 + 1440 ^ ^ + 810 tt^ ^ x + 405 tt^ x^ + 720 x 

- 1440 - 896 TT^ . 

Eqs. (p])-(|52|) provide the full metric up to and including the order 0{fi^). 

Making use of the metric functions, one can furthermore find the series ex- 
pansion for other quantities of interest: 



V2 f 8 1 
~ V 9^ ^ 16 



0(^9/2) (54) 



and 



/ 163^2 77^ 7V2?73 3V2r] ^ y/2r] ^5/2 



„5 



V 4807r2 487r2 327r2 

9V2?y 163V2r]-' iV^rj^ 59?/^ %/2 1117ry^%/2 

" 1287r2 ^ 19207r2 1287r2 36^ ^ 134407r2 

777372 \/2 7y'5 25\/2 7/' 



127r4 47r4 SStt^ 



0(m^/2). (55) 



Related to the fact that the first four /„ can be calculated without making 
use of the Jacobi inversion problem, Bardeen and Wagoner were able to de- 
termine them in Q. Taking into account that the relativistic parameter used 
in said paper is somewhat different, |^ the coefficients calculated by Bardeen 
and Wagoner were compared with eqs. (|3^)-(|3^) and found to be in agreement. 
Using the knowledge that these first coefficients are correct, a further twelve /„ 
have been proved correct by showing that they satisfy the Ernst equation up 
to the relevant order. The further eight coeffecients that were explicitly calcu- 
lated could not undergo the same test due to computer limitations, but were 
calculated using the same iteration scheme, whence their correctness is all but 
certain. For more details see Q. 

While a direct comparison of the PN expansion of the Ernst potential to 
the numerical results of Bardeen and Wagoner is not feasible for an arbitrary 
point in space, various quantities in the disc can be compared quite readily. 
In particular, the analytical version of Table 2 in § was drawn up, presenting 

^Bardeen and Wagoner chose to use the relativistic parameter 7 defined by 7 = 1 — e^". 
The series forms of the functions 7 = 7(^1) and fi = fJ,{'y) can be calculated using the furK^tion 
Vb = Vb(/i) (see H) and are listed up to the orders 0{^^^) and O^-y^^) respectively in |H|. 
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the series for the square of the linear velocity, u^, for some particle of the disc 
as measured in the locally nonrotating frame of reference. It turns out that 
the numerical results, which are listed to the 5*^^ decimal place, are perfectly 
correct up to the order 0{-f^). ^ As to be expected, the expansion coefficients 
differ more and more from the analytic values as one moves to higher orders 
of the series. The last entry in Table 2 of Q deviates from the analytic result 
by approximately 25%. In all however, the numerical results in Bardeen and 
Wagoner are highly accurate and seldom differ from the analytic ones by more 
than a fraction of a percent. 

4.2 Convergence 

As could have been expected right from the outset, the convergence of the PN 
approximation depends to a great extent on position, i.e. on the coordinates ^ 
and 77. For example, on the rim of the disc itself, the Ernst potential is given by 
the analytic expression / = 1 — /i/2, whereas in general, the potential is given by 
an infinite series in ^. Plots similar to that of Fig. || strongly indicate that the 
series converges for spatial infinity ^ 00) and has a radius of convergence of 
precisely /io in this limit. Evaluations of the series at various (arbitrary) points 
in space indicate, moreover, that it converges everywhere. The series tends 
to converge quite quickly for small \r]\, but rather erratically for \ri\ nearing 1 
(77 = ±1 corresponds to the axis of symmetry). 

The Pade approximant, which represents a truncated series via a rational 
polynomial expression, proves to be a highly advantageous alternative to the 
series, in particular if the original series converges slowly, erratically or not at 
all. Using this approximant, the erratic convergence of the PN expansion for 
the Ernst potential near the axis or the particularly slow convergence of the 
dimensionless quantity Qpo can be largely circumvented. Frequent use of the 
Pade approximant was made in the work of Bardeen and Wagoner Q . In this 
paper however, sparing use of it is made so as better to concentrate on the PN 
expansion itself, but a comparison of the Pade approximant to the standard PN 
series can be found in Table |^. Let it be mentioned that the proper surface 
mass density (Fig. ^ and the ergospheres (Figs. ^ and |^) can be found to 
extremely high accuracy using the Pade approximant. A general discussion of 
the Pade approximant can be found in and a consideration of its use in the 
PN expansion in Q. 

In Figure 0, the ratio of consecutive coefficients from the real and imaginary 
parts of the Ernst potential has been plotted to illustrate the convergence of the 
series for a chosen point in space. In order to have a measure for the leading 
terms of the expansion, we took /o to be 1 and /_i to be i. Note that the 
series is guaranteed to converge at a given value of /i, say fi — iic, so long as the 
condition lim„^oo \ fn+i/ fn-i\ < l//^c is met. Thus the series always converges 
in the Newtonian limit fi and converges in the limit /i — > /io so long as 
lim„^oo |/„+i//„-i| < l//io ~ 0.21 holds. 

^The results concerning v'^ were calculated using the parameter 7 in order to facilitate the 
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Figure 1: The ratio of consecutive real and imaginary coefficients from the 
expansion of / taken at a point ^ = 0, rj = 3/5 in the disc. The sohd hne 
represents the real coefficients and the dotted line the imaginary coefficients. 



In Fig. |l| one sees that the series converges quite quickly for this point on 
the disc, that /i plays a fairly important role, but that /a and higher coeffi- 
cients contribute little to the potential. The fact that the ratio of consecutive 
coefficients of the series is negative means that these coefficients alternate sign, 
so that one always has an upper and lower bound on the value for (the real 
and imaginary parts of) /. Note that not all points in space exhibit such clear 
convergence as that of Fig. [|. 

Not only does the applicability of the series depend upon the point in space 
that is chosen, but of course upon the value of ^, as well. It goes without saying 
that the PN approximation reflects the full field equations fairly accurately for 
small values of /i, but we shall see further, as did Bardeen and Wagoner much 
to their surprise, that the approximation can be utilized to great benefit even in 
many highly relativistic situations. Looking at Table |l|, one can see that the PN 
approximation of order ©(/i^^/^) returns values for / that are correct to at least 
9 decimal places at /i = 1/2 (the numerically evaluated potential, determined 
using the analytic solution, is correct at least as far as it is given in Table |l|). 
The highly relativistic case /x = 3 can be handled more than satisfactorily by 
the PN approximation, whereby the potential on the axis is valid to only four 
decimal places. The accuracy of the series approximation is astonishingly high 
even for the limit fj, fj.^. Although the results on the axis are reasonably 

comparison to the work of Bardeen and Wagoner. 

^It should be noted that po tends to zero in this hmit. The extreme Kerr metric is obtained 
for + C'^ 7^ 0, whereas the analytic solution shows that finite values for ^, as considered 
here, lead to a different (not asymptotically flat) spacetime (see [BI). The PN approximation, 
on the other hand, erroneously yields an asymptotically flat spacetime in this limit, at least 
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M = 1/2 


^ = 3 


= Mo 


V = 1/2 


0.772687415 - 
0.022438224i 
0.772687415- 
0.0224382244 


-0.281941980 - 
0. 3528103991 
-0.281941981 
0.352810399i 


-0.861419387 - 
0.6849241174 
-0.861419898- 
0.6849243964 


On the axis 
e = l, 
'7=1 


0.820284295 - 
0.032540154i 
0.820284295- 
0.0325401544 


0.127426794 - 
0.560328316i 
0.127409032- 
0.5603159544 


0- 

4 

0.006959726 
0.9977838114 


In the disc 

e = o, 

7? = 3/5 


0.675204587 - 
0.041948129i 
0.675204587- 
0.0419481294 


-0.626305938 - 
0.409718317?; 
-0.626305488- 
0.4097180904 


-1.304973284 - 
0.6167998884 
-1.304859086- 
0.6167285384 



Table 1: The Ernst potential / according to a numerical evaluation of the exact 
solution in comparison with the results of the PN approximation (bold-faced) 
up to the order ©(^^^Z^) for various points in space and three values of fi. 

good in this limit (accurate to within 1%), one can obtain much more accurate 
results using the Pade approximant. Table ^ provides a comparison of the value 
of / for a point along the axis returned by the PN series to that of the "diagonal" 
Pade approximant, whereby diagonal is meant to indicate the approximant for 
which the polynomial in the numerator is of the same order as the polynomial 
in the denominator. The values in this table should be compared to the "exact" 
value found in the cell of the 2"'^ column and 2"^^ row of Table |l|. It turns out 
that the values returned by the Pade approximant are so accurate as to provide 
a viable alternative to numerical methods for all practical purposes. ^ Using 
Table |^ one can see, moreover, that the nature of the convergence of the Pade 
approximation is more uniform and hence more predictable than that of the 
PN approximation. Thus one could have come up with a good estimate for 
the accuracy of the Pade values in the table, but not for the PN values. The 
PN value of order 0{fi^^^^), for example, is scarcely more accurate than that 
of order 0(/i^'^/^) whereas the corresponding Pade values have won a further 8 
decimal places of accuracy. As a rule of thumb, it appears one can rely on the 
Pade approximant to gain one decimal place of accuracy for each increase of y/Ji 
in the order of the polynomial. 

A pictorial impression of the convergence of the PN series can be gleaned from 
Fig. |. The proper surface surface mass density dp, which can be calculated 
from ap = ae^^^ where a is given by eq. (23) of ||], is divided by Vt to give a 
dimensionless quantity and plotted as a function of the normalized radius p/po- 
A discussion of the curve itself, which lies outside the scope of this paper, can 
be found in Q . What is of interest here is the way in which the curves approach 
the curve (f). Although the tendency to converge can be seen quite clearly, the 
discrepancy between curves (e) and (f) is fairly large. This is primarily due to 

so long as the series is finite. 

* However, as /.i — > /.io, one has to be careful for large 
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0(^^21/2) 


0(^25/2) 


PN 


0.150938555- 
0.564116343i 


0.127710083- 
0.556936459i 


0.126549231- 
0.561010096* 


0.127627187- 
0.5602818771 


0.127409032- 
0.560315954j 


Fade 


0.140716023- 
0. 5624266874 


0.127381473- 
0.560314135i 


0.127426866- 
0.560328310j 


0.127426794- 
0.560328316i 


0.127426794- 
0.560328316j 



Table 2: The Ernst potential / for /i = 3 at the point ^ = 1, ?; = 1 according 
to the PN approximation and the Padc approximant for different orders of the 
series (cf. 2"^^ row, 2"^^ column of the previous table). 

the slow convergence of the dimensionless quantity flpo. For example, the PN 
approximation of order 0{^'^^^'^) returns a value of approximately 0.220 instead 
of 0.213 for fipo at /i = 3 and aX fi — 4.5 the PN approximation yields the grossly 
erroneous value of 1.46 as compared with the correct value, 1.01 x 10^^. With 
the Pade approximant, this problem vanishes and the curves for the proper 
surface mass density are indistinguishable from the exact curves even at the 
order ©(/.t^/^). 

4.3 Ergospheres 

The results of Table |^, i.e. the fact that the PN approximation returns a highly 
accurate value for / even in very relativistic situations, justify the use of this 
approximation even for the consideration of purely relativistic phenomena. One 
such phenomenon, the ergosphere, dramatically highlights differences between 
the Newtonian and Einsteinian theory. An ergosphere is a region in which the 
Killing vector characterising the stationary nature of the spacetime becomes 
spacelike (the metric function e^^ is negative in this region) so that nothing can 
remain still relative to an observer at infinity — everything is dragged along by 
the rotation of the disc. 

Figures ^ and |4| depict the ergosphere for two values of /i and calculated for 
three different orders of the PN approximation (please see |]lOj for a comparison 
with ergospheres found by evaluating the analytic solution numerically). One 
can see in Fig. || that the ergosphere for /i = 2, as given by the 6*^ PN approx- 
imation, is nigh on correct, so that higher order terms serve only to refine fine 
details. In the case of = 3, however, (Fig. ^) one can see marked improvement 
as one moves toward higher orders of the approximation. 

As ^ increases, the ergosphere moves closer and closer to the axis of rotation 
where the PN approximation is no longer reliable. This, compounded with the 
fact that the increasingly poor convergence of the dimensionless quantity fipo as 
/i nears /iq results in scaling problems, leads one to assume that the PN approx- 
imation is ill-suited for rendering the ergoshpere for values of /i approaching /iq- 
Due to the extreme accuracy of the Pade approximant however, one can indeed 
divine the nature of the ergospheres even in this limit. 
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Figure 2: The radial distribution of the proper surface mass density for ^ = 3 at 
different orders of the PN approximation. The curves were created by employing 
the approximation up to the order (a) ©(/z^/^), (b) 0{ji^^/'^), (c) ©(/i^'^/^), 
(d) ©(/i^^/^) and (e) 0{ii^^^'^). The curve (f) was created by evaluating the 
exact solution numerically to extremely high accuracy. 
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Figure 3: The ergospheres as calculated using the 6**^, 8**^ and 12*'' PN ap- 
proximation for /i = 2. The outline represents the curve along which e^^ is zero 
and should be thought of rotated about the axis of symmetry. The inside of the 
resulting torus-like figure of revolution is then the ergosphere. 
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Figure 4: The ergospheres as calculated using the 6* , 8* and 12 PN ap- 
proximation for /i = 3. The outline represents the curve along which e^^ is zero 
and should be thought of rotated about the axis of symmetry. The inside of the 
resulting torus-like figure of revolution is then the ergosphere. 



5 Conclusion 

It was possible, due to the existence of an analytic, global solution for the 

axially symmetric, stationary, rigidly rotating disc of dust, to come up with an 
iteration scheme to calculate an arbitrary coefficient in the PN approximation 
of this solution. This work amounts to the analytic analogue of numerical work 
published by Bardeen and Wagoner in 1971. The explicit calculation of the 
expansion coefficients of the Ernst potential was carried out to the 12"^ PN 
level (i.e. It turns out that the numerical results from Bardeen and 

Wagoner are highly accurate for lower orders of the expansion and remain quite 
good even for higher orders. 

The supposition that the PN approximation could be used for many highly 
relativistic situations was confirmed for the disc of dust, and the very accurate 
renderings of the ergospheres obtained using the PN approximation attest to 
the fact that physically meaningful, relativistic phenomena can be studied in 
some cases using tools tailored to other applications (i.e. to the consideration of 
the Newtonian regime) . It turns out that the position in space to be considered 
contributes significantly to the convergence of the PN approximation, but not 
to the convergence of the corresponding Padc approximant. It was found that 
the PN approximation is unreliable for large ^ in the neighbourhood of the axis 
of rotation even when the approximation is quite accurate elsewhere for the 
same value of /i, and it can be made highly accurate for all points in space at 
all values of /i by applying the Pade approximation methods. If one is careful 
to take into account its limitations, one can use the PN approximation to great 
advantage for quick and accurate calculations and is bolstered in the opinion 
that approximation methods can sometimes be extended to applications outside 
their a priori guaranteed region of validity. 
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